
% inter-arrival times (PH distribution)
alpha = [1/2 1/2];
T = [-1 0;
     0  -2];

% service times (PH distribution)
beta = [1 0];
S = [-10 10; 
     0  -10];

[ql,wait_alpha,wait_T]=Q_CT_PH_PH_1(alpha,T,beta,S);